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Abstract 

In  this  study,  both  solid-phase  and  fluid-phase  temperatures  inside  a  porous  cathode  of  a  polymer  electrolyte  fuel  cell  in  contact  with  an 
interdigitated  gas  distributor  are  predicted  numerically.  The  porous  cathode  consists  of  a  catalyst  layer  and  a  diffusion  layer.  The  heat  transfer  in 
the  catalyst  layer  is  coupled  with  species  transports  via  a  macroscopic  electrochemical  model.  On  the  other  hand,  in  the  diffusion  layer,  the  energy 
equations  based  on  the  local  thermal  non-equilibrium  (LTNE)  are  derived  to  resolve  the  temperature  difference  between  the  solid  phase  and  the 
fluid  phase.  As  for  the  species  transports,  the  Bruggemann  model  is  employed  to  describe  the  effective  diffusivities  of  the  oxygen  and  water  vapor 
in  the  porous  cathode.  Results  show  that  the  wall  temperature  decreases  with  increasing  the  intrinsic  heat  transfer  coefficient.  As  the  intrinsic  heat 
transfer  coefficients  increase,  the  porous  electrode  becomes  local  thermal  equilibrium  with  a  strong  thermal  interaction  (heat  transfer)  between  the 
solid  and  fluid  phases.  Under  the  conditions  of  high  intrinsic  heat  transfer  coefficients,  the  temperature  difference  between  the  solid  matrices  and 
the  reactant  fluids  are  negligible. 

©  2006  Elsevier  B.Y.  All  rights  reserved. 
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1.  Introduction 

Thermal  transports  in  a  polymer  electrolyte  fuel  cell  (PEFC) 
have  received  a  particular  attention  due  to  their  strong  effect 
on  the  fuel  cell  durability.  The  electrode  of  a  PEFC  is  a  typi¬ 
cal  porous  medium.  In  this  medium,  the  heat  dissipation  by  the 
electrochemical  reaction  in  the  active  layer  of  a  fuel  cell  should 
be  removed  thoroughly  to  sustain  proper  operation. 

There  are  two  tasks  in  the  thermal  management  of  a  low- 
temperature  fuel  cell.  The  first  is  to  optimize  the  operating 
temperature  range  that  avoids  the  local  hot  spot  as  well  as  the 
water  flooding.  The  second  task  is  to  keep  a  uniform  temper¬ 
ature  distribution  that  avoids  the  thermal  stress  in  the  porous 
electrode.  Therefore,  a  superior  thermal  management  becomes 
critical  for  safely  operating  and  prolonging  a  fuel  cell. 

However,  most  investigations  about  the  porous  electrode  of  a 
fuel  cell  did  not  consider  thermal  transports  [1-5]  and  only  a  few 
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discussed  the  thermal  behaviours  in  the  porous  electrode  based 
on  the  local  thermal  equilibrium  (LTE)  assumption  [6-9].  The 
LTE  approach  is  simple  and  straightforward  but  is  only  appli¬ 
cable  if  the  temperature  difference  between  the  solid  and  the 
fluid  phases  is  very  small.  The  so-called  one-equation  model 
is  employed  to  resolve  the  phase-mixing  temperatures  in  the 
porous  medium.  When  the  temperature  difference  is  signifi¬ 
cant,  it  becomes  a  local  thermal  non-equilibrium  (LTNE).  The 
one-equation  model  is  invalid,  and  should  be  modified  with  the 
two-equation  model,  which  treats  the  solid-phase  temperature 
and  the  fluid-phase  temperature  separately.  An  analysis  of  heat 
transfer  in  a  porous  medium  based  on  the  two-equation  model  is 
more  involved  because  it  requires  interstitial  information  inside 
the  porous  medium  such  as  heat  transfer  coefficient  between  the 
fluid  and  solid  phases  and  the  interfacial  surface  area  [10-12]. 
Due  to  this  difficulty,  most  investigations  about  the  porous  elec¬ 
trode  of  a  fuel  cell  resorted  to  using  the  one-equation  model  to 
obtain  the  phase-mixing  temperatures. 

Studies  discussing  the  phase-temperature  difference  between 
the  reactant  fluids  and  the  solid  matrices  inside  the  porous  elec¬ 
trode  of  a  polymer  electrolyte  fuel  cell  are  rather  sparse  [13,14], 
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Nomenclature 

Bi 

Biot  number 

Cu20 

water  vapor  mole  concentration  (molm-3) 

CH20,ref 

water  vapor  mole  concentration  at  inlet  (mol  m-3 ) 

Cr2o 

normalized  water  vapor  concentration,  ch2o /Got 

^H20,ref 

normalized  water  vapor  mole  concentration  at 
inlet,  Cn20,ref/ <4ot 

co2 

oxygen  concentration  (molm-3) 

c02,ref 

oxygen  mole  concentration  at  inlet  (molm-3) 

Co2 

normalized  oxygen  concentration,  co2/ctot 

f-02,ref 

normalized  oxygen  mole  concentration  at  inlet, 

c02,ref/  Got 

CP 

specific  heat  at  constant  temperature  (J  kg  - 1  K- 1 ) 

ctot 

total  mole  concentration  of  the  reacting  fluid 
(molm-3) 

d 

pore  diameter  of  the  porous  medium  (m) 

£>h2o 

binary  diffusivity  of  water  vapor  in  the  oxygen 
(m2  s-1) 

^H20,eff 

effective  diffusivity  of  water  vapor  in  the  gas  dif¬ 
fusion  layer  (m2  s-1) 

Do2 

binary  diffusivity  of  oxygen  in  the  water  vapor 
(m2  s-1) 

D(h.eif 

effective  diffusivity  of  oxygen  in  the  gas  diffusion 
layer  (m2  s-1) 

Da 

Darcy  number 

Ec\ ,  Ec2 

coefficients  in  Eq.  (23) 

E 

Faraday’s  constant  (96487  Cmol-1) 

hy 

interstitial  heat  transfer  coefficient  in  the  porous 
medium  (W  m  3  K-1) 

I 

current  (A) 

k 

thermal  conductivity  (W  m-1  K-1) 

L 

length  of  the  computational  domain,  160  [xm  (m) 

m 

mass  flux  (kgm-2  s-1) 

P 

pressure  (Pa) 

P 

non-dimensional  pressure 

Pr 

Prandtl  number 

Q 

heat  dissipation  by  the  electrochemical  reaction 
(W) 

Re 

Reynolds  number 

7?mi  ?  7?m2  coefficients  in  Eq.  (24) 

Rni,  7?n2  coefficients  in  Eq.  (25) 

Rq 

ohmic  resistance  (£2) 

Sy 

surface  area-to- volume  ratio  (m-1) 

Sc 

Schmidt  number 

St 

Stanton  number 

T 

temperature  (K) 

U,  V 

velocity  components  in  the  x,  and  y  directions, 
respectively  (ms-1) 

u,v 

non-dimensional  velocity  components  in  the  x, 
and  y  directions 

W 

channel  width,  40  |xm 

ws 

shoulder  width,  80  p,m 

x,y 

coordinate  system  (m) 

X,Y 

non-dimensional  x  and  y  coordinate 

Greek  symbols 

a  i ,  c*2  coefficients  in  Eq.  (11) 

8  thickness  of  the  diffusion  layer,  40  jxm  (m) 

8C  thickness  of  the  catalyst  layer,  10  jxm  (m) 
s  porosity  of  the  diffusion  layer 

£c  porosity  of  the  catalyst  layer 

0  cathodic  overpotential  (V) 

6  non-dimensional  temperature 

k  permeability  (m2) 

p  density  (kg  m- 3 ) 

r  tortuosity 


and  this  has  motivated  the  present  research  to  develop  a  two- 
equation  model  to  predict  the  fluid  and  solid  temperatures  inside 
the  porous  electrode  of  a  polymer  electrolyte  fuel  cell.  The 
effects  of  the  intrinsic  heat  transfer  coefficient  of  the  porous  elec¬ 
trode  on  the  heat/mass  transfer  mechanisms  are  examined  and,  as 
a  result,  thermal  paths  in  a  porous  cathode  of  polymer  electrolyte 
fuel  cell  are  clearly  identified.  It  would  be  of  assistance  for  fur¬ 
ther  accurate  analyses  of  the  fuel-cell  thermal  performance  by 
considering  the  thermally  dependent  physical  properties  inside 
a  polymer  electrolyte  fuel  cell. 

2.  Numerical  model 

Fig.  1  (a)  shows  a  sectional  view  of  a  porous  cathode  in  contact 
with  an  interdigitated  flow  field  of  a  polymer  electrolyte  fuel 
cell.  The  present  computational  module  is  limited  in  the  module 
shown  in  Fig.  1(b).  The  oxygen  reduction  reaction  (ORR)  on  the 
cathode  is  considered  in  the  present  model. 

02(g)  +  4H+  +  4e“  **  2H20(g)  ( 1 ) 

The  oxygen-water  vapour  mixture  enters  the  porous  cathode 
form  the  module  inlet  (inlet  channel),  and  then  trans verses  the 
diffusion  layer  to  the  catalyst  layer.  The  electrochemical  reaction 
occurring  in  the  catalyst  layer  consumes  the  oxygen  (Co2)  and, 
meanwhile,  produces  the  water  vapour  (Ch2o)«  In  addition,  the 
heat  ( Q )  is  generated  by  the  overpotential  and  irreversibility 
during  the  electrochemical  reaction.  Either  the  fluid  or  the  solid 
removes  it.  The  assumptions  in  this  model  include:  (i)  mixtures 
are  ideal  gas;  (ii)  flow  is  steady,  laminar,  and  incompressible;  (iii) 
thermal  physical  properties  are  constant;  (iv)  porous  electrode  is 
homogeneous  and  isotropic;  (v)  water  exits  as  vapour  only  and 
(vi)  both  inlet-fluid  temperature  and  rib-surface  temperature  are 
uniform. 


2.7.  Governing  equations 


In  the  porous  layers,  the  steady  volume-averaged  continuity 
and  momentum  equations  that  govern  fluid  flow  are  given  by 


d(pfu)  d(pfv) 


dx 


dy 
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Fig.  1.  (a)  Schematic  drawing  of  porous  electrode  of  the  interdigitated  flow  field,  (b)  Configuration  of  the  computational  domain. 


d2u 

d2u N 

J  + 

(3) 

dx2 

1  dy 2y 

d2v 

d2v\ 

+  Sy 

(4) 

dx2 

1  dy2) 

du  du  dp 

PfuC - 1"  PfvC~  =  ~  X - l~  M 

ox  dy  ox 

dv  dv  dp 

Pfu  ~ — I"  PfF—  =  —  - — ^  P 
dx  dy  dy 

The  source  terms  in  Eqs.  (3)  and  (4)  relate  to  the  traditional 
Darcy  dissipation.  It  is  used  to  account  for  the  pressure  drop 
due  to  viscous  friction  at  the  surfaces  of  carbon  fibbers,  i.e., 
Su  =  —spu/K  and  Sy  =  — spv/K ,  respectively.  As  for  the  energy 
equations,  the  two-equation  model  is  used  to  describe  the  ther¬ 
mal  behaviour  in  the  diffusion  layer  [12].  The  energy  equations 
for  solid  and  fluid  phases,  respectively,  are: 


,  dZTs  dzTs 

0  =  *s,eff  (  —2  + 


dTf 

( PCp\u —  +  (pcp\v 


-  hy(Ts  -  7f)  - 


(5) 


d2rf 


+ 


dx2  dy2 
T  hv(Ts  —  7f) 


(6) 


The  last  two  terms  of  Eq.  (5)  mean  the  convective  heat 
delivered  to  the  reactant  fluids  from  the  carbon  fibers  and  the 
Joule-heating  effect,  respectively.  The  last  term  of  Eq.  (6)  rep¬ 
resents  the  heat  received  by  the  reactant  fluids  from  the  carbon 
fibers. 

The  effective  thermal  conductivities  for  the  solid  phase  and 
the  fluid  phase  are  respectively  defined  as 


^s,eff  —  (1 
£f,eff  =  Okf 


(7) 

(8) 


In  the  catalyst  layer,  physically,  the  electrochemical  reaction 
occurs  at  the  interface  of  reactant  (fluid)  and  catalyst  (solid) 
under  a  fixed  temperature.  Therefore,  the  fluid  and  solid  phases 
in  the  catalyst  layer  have  the  same  local  temperatures,  i.e., 


dTf  dTf  f 

(Pcp) ^  (PCP^ f^^T  —  ^c,eff  ( 


3  2T{  327f 
dx2  +  dy2 


+  JT  3 


(10) 


The  source  term  jjr]  represents  the  overpotential  heating  by  the 
electrochemical  reaction.  According  to  the  Butler- Volmer  cor¬ 
relation  [15],  the  relationship  between  the  local  transfer  current 
density  and  the  reactant  concentrations  can  be  described  as  the 
following  equation. 

2 


J  T  =  Oil 


c  02 


(  cu2o 

a2[ 


(11) 


C02  ,ref  J  \  OH2  O  ,ref 

a  i  and  d2  are  electrochemical  coefficients  depending  on  the 
exchange  current  density,  the  overpotential  on  the  electrode  sur¬ 
faces,  and  the  reacting  temperature.  It  is  known  that  the  kinetics 
of  oxygen  reduction  reaction  is  related  to  the  temperature.  Since 
the  temperature  variation  in  reacting  layer  for  the  PEFC  is  small, 
the  temperature  effect  on  the  kinetics  of  ORR  is  negligible.  Con¬ 
sequently,  the  coefficients  a\  and  d2  are  regarded  as  constants 
in  the  present  model. 

The  species  transports  of  the  oxygen  and  water  vapor  in  the 
porous  cathode  can  be  governed  by  the  following  equations: 

dx2 


u 


dco2 dc02 


U 


dx 

dcii20 

dx 


+  v 


+  v 


dy 

dcii2o 


,  dzco0  d2co0 
A)2,eff  (  ~~X~S^L  + 


dy2 


+  5o2 


dy 


7JH2Q,eff 


d2CH2o\  , 

J 


In  the  catalyst  layer,  the  sources  terms  for  oxygen  trans¬ 
port  and  water- vapor  transport  are  So2  =  Jt/4E  and  Sh2o  = 
—j^/2F,  respectively.  In  the  diffusion  layer,  they  are  0.  The 
effective  diffusivities  of  the  oxygen  and  water  vapour  in  the 
porous  electrode,  Do2,eff  2ind  Dn?o,eff?  follows  the  Bruggemann 
model  [16],  i.e., 


(9)  Ao2,eff  =  sT  Dq: 


Tf  =  Ts 


(14) 


Table  1 

Governing  equations  of  the  present  model 
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Equations 


Sc J> 


Diffusion  layer  Catalyst  layer 


Energy,  diffusion  layer 


Energy,  catalyst  layer 


Species 


Mass 

dU  , 
dx  “c 

dV 

dY 

=  s& 

3  U 

3  U 

3  P 

1 

Momentum 

U 

+  V 

— 

— 

+ 

dX 

37 

dX 

Re 

dV 

3V 

3  P 

1 

u 

+  V 

—  . 

+ 

dX 

37 

97 

Re 

d2U  d2U 
+ 


dX 2 

d2v 


+ 


dY2 

d2V 


dX2  dY 2 


+  S& 


+  s@ 


Non-dimensional  parameters 


0  = 


d2e. 


+ 


d2es 


dX2  3 Y2 


+  s& 

1 


86f  d6f 

U—  +  V—  ■ 

dX  dY  RePrt  ff 


0S  =0f 


1 


dOf  dOf 

U  —  +  7  — - 
dX  dY  RePreff 


3Co,  dC  o, 

U—fA  +  v  2 


1 


dX 


dY  ReSco, 


320c  32@f 

iF+  W  1  +  s* 


92ef  d2ef 
ax2  +  W]+S0 

d2C0l  32C0, 


+  s@ 


jj  9Ch2o  y  9Ch2o 


1 


dX2  dY2 

92Gh2o  ,  92Gh2o 


dX 


x 

X  =  . 

8 


dY 


Y  =  -, 
8 


ReScu^o  \  dX2 


+ 


dY2 


+  S$ 


m 


Co2,ref  — 


c02,ref 


Pr eff  - 
^M1  = 
^N2  - 


Ctot 
/x(Cp)f 

aq32 

4  E ( Ctot  Do2  ,eff 

aid2 

2FctotDH20,eff  ’ 


Md  = 

Gh2o  = 


Pf  Ac  ’ 

ch2o 

Ctot 
2 


Bi  = 


hsSvd 

ks,eff 

^M2  - 


U 

U=  — , 
wd 

Ch20,  ref  — 
hsSyd 


V 


V=—, 

u& 


CH20,  ref 


Ctot 


0 


1 


ReDa 
1 


■eU 


■eV 
ReDa 

Bm  -  Os) 

sm s  -  of) 


0 


0 


0  = 


p 

p  =  — — 
Pf^d 
T  -Tr 


7f,in  -  Tr  ’ 


= 


(pCp)fMd  ’ 


Scq2  — 


v 


S( J> 


1 


ReDa 

1 

ReDa 


sCU 


sCV 


1 


RePre  ff 

1 

ReSco2 

1 


RcScH20 

C02 


[EciCqo  —  £c2(CH9o)  ] 


[^MlCoo  —  /?M2(Ch9o)  ] 


[RmCo2  —  Rm(C^o)  ] 


co?  = 


Re  = 


Ctot 

Udd 


&2  d" 


K 


Rk  = 


4  F Ctot  Dq2  eff 
&f,eff 


^s,eff 


Ec\  = 


Da  ^ , 
a  i  r]dz 


D02  ,eff 

/?N1  = 


^s,eff(^f,in  Tr)Co2,ref 


,  5ch2o  - 

ai^2 

2Ectot  Dn2o,eff  ’ 
EC2  = 


V 

V 


Du20,  eff 


0L2dd" 


&s,eff(Tf,in  ^r)GH2Oref 


Table  2 

Porous  electrode  data  sheet 


Property  name 

Physical  propertiesa 

Material 

Carbon  fiber 

Bulk  density,  p 

1100  kg  m-3 

Specific  resistance,  tr-1 

115  [x£2m 

Thermal  conductivity,  k 

1.71  W  m-1  K-lb 

Geometric  properties 

Diffusion  layer  porosity,  e 

48% 

Catalyst  layer  porosity,  £q 

42% 

Diffusion  layer  thickness,  d 

400  pan 

Catalyst  layer  thickness,  3c 

100  [Jim 

Diffusion  layer  tortuosity,  z 

1.5 

Catalyst  layer  tortuosity,  z 

1.5 

Average  pore  diameter,  d 

33  [Jim 

Diffusion  layer  permeability,  kc 

1.57  x  10"12m2 

Catalyst  layer  permeability,  kqc 

1.02  x  10“12m2 

Flow  properties 

Interstitial  heat  transfer  coefficient,  hy 

1.0  x  103  to  1.0  x  109  W  m~3  (assumed) 

Surface-to-volume  ratio,  Sy 

1.0  x  103  m-1  (assumed) 

Pore  velocity,  Ud  for  /=  2.0  A  at  stoichiometric  flow  ratio  5 

0.478  ms-1 

a  Data  from  manufacture  (Toray). 

b  Oxygen  at  25  °C  with  15%  RH;  p  =  0.15  atm;  plate  thickness,  2.5  cm;  51  Lm~2  s_1. 
c  ASTM  E1530. 
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£>H2o,eff  =  £t£>h2o  (15) 

The  above  governing  equations  can  be  deduced  to  dimen¬ 
sionless  forms  of  using  proper  non-dimensional  parameters  in 
Table  1.  Note  that  according  to  the  results  obtained  by  using 
the  order-of-magnitude  analysis  [17],  the  Stanton  number  of  a 
typical  diffusion  layer  is  about  1 . 

2.2.  Boundary  conditions 

The  boundary  conditions  for  of  the  present  computational 
domain  are  summarized  as  follows: 

(1)  At  the  module  inlet, 

p  =  Pin ,  Of  =  0m,  ^  =  Bi(Os  -  Of), 

OX 

Co2  =  0.9,  CH2o=0.1  (16) 

The  relative  humidity  of  the  inlet  gases  is  about  50%  RH  at 
the  operation  temperature  of  60  °C. 


Fluid-Phase  Temperature  Max:  1.183 


Min:  1.00 


Velocity  Vectors 


1.6 

1.5 

1.4 

1.3 

1.2 

1.1 

1 

0.9 

0.8 

0.7 

0.6 

0.5 

0.4 

0.3 

0.2 

0.1 

0 


iff 


t  T  T 


\  t  1 


T  1  \ 


t  \  ^ 


^  ^  h*  1  jt — ^  ^ ^ — y* — ^ 


t  // 

ntmA 

&■ 


'VVVVvVVV'WV' 
'VvVWVY  VWV' 

d/ 


<  < 


-0.2  -0.1  0  0.1  0.2  0.3  0.4  0.5  0.6 


Fig.  3.  Comparison  of  fluid-phase  and  solid-phase  temperature  distributions 
inside  the  porous  electrode  for  Re  =  6,  Da  =  9.8  x  10“6,  and  St  =  0.74. 


Fig.  2.  Velocity  distributions  inside  of  porous  cathode  of  a  polymer  electrolyte  Fig-  4.  Comparison  of  fluid-phase  and  solid-phase  temperature  distributions 
fuel  cell.  inside  the  porous  electrode  for  Re  =  6,  Da  =  9.8  x  10’6,  and  St=  14.73. 
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(2)  At  the  rib  (current  collector)  surfaces, 

dP 


dX 

dof 

dx 


=  o, 


=  sm  -  of), 


U  =  V  =  0, 
9Cq2 


dx 


os  =  i, 

dCff2o 

'  dx 


=  o 


(17) 


(3)  At  the  interface  of  the  membrane  and  the  catalyst  layer 

(*=-0.1), 


U  =  V  =  0, 

dP  dos  dof 


dCo2  3Ch2o 


dx  dx  dx  dx 

(4)  At  the  symmetric  planes, 

dP  du  dv  dos  dOf 


dX 


=  0 


(18) 


dC  o2  9Ch2o 


dr  dr  dr  dr 


dr 


dr 


dr 


o 

(19) 


(5)  At  the  module  outlet 

dof 


p  =  p 


out' 


dx 


=  o, 


90s 

dx 


=  Bi(0s  -  Of), 


dC  o2  9Ch2o 


dX 


dX 


=  0 


(20) 


2.3.  Numerical  approaches 


The  governing  equations  are  numerically  solved  by  the  finite- 
element  method.  It  uses  the  Broyden’s  method  with  an  LU- 
decomposition  pre-conditioner  to  solve  the  nonlinear  equations 
iteratively.  A  penalty  term  is  employed  for  pressure  to  reduce 
continuity  errors.  Thus,  there  is  a  continuous  part  of  the  pressure 
and  piecewise  constant  part  providing  and  extra  DOF  (degree  of 
freedom)  for  pressure  on  each  element.  It  uses  Newton-Raphson 
iteration  to  solve  the  close-coupled  groups  (velocity,  pressure, 
temperature,  concentration  and  electricity)  and  uses  the  frontal 
algorithm  (Gaussian  elimination)  to  solve  the  linearized  system 
of  equations  for  each  iteration.  Variable  grids  in  the  v-direction 
and  uniform  grids  in  the  y-direction  were  employed  in  the  present 
study.  Solutions  are  considered  to  be  converged  when  all  residual 
sources  (including  mass,  momentum,  temperature,  species  and 
charge)  were  less  than  1 .0  x  10~6.  A  typical  simulation  requires 
about  60  min  of  central  processing  unit  time  on  a  Pentium  IV 
2.80  GHz  PC. 


3.  Results  and  discussion 

3.1.  Velocity  distributions 

Fig.  2  shows  the  distribution  of  flow  velocity  vectors  in  the 
porous  electrode.  The  Reynolds  number  and  Darcy  number  are 
fixed  at  Re  =  12  and  Da  =  9.8  x  10-6,  respectively.  It  is  seen  that 
the  velocities  are  higher  in  the  diffusion  layer  than  in  the  catalyst 
layer.  This  is  because  the  diffusion  layer  accompanies  a  higher 
porosity  and  a  higher  permeability.  In  addition,  the  velocities 
near  the  shoulder  surface  are  higher  as  compared  to  those  adja¬ 
cent  to  the  catalyst  layer.  This  is  reasonable  because  the  flow 


paths  from  the  inlet  channel  to  the  outlet  channel  are  shorter  for 
the  flow  closer  to  the  shoulder  surface.  A  shorter  flow  path  has 
a  lower  pressure-drop  penalty,  which  allows  more  fluids  to  pass 
through  (Table  2). 

3.2.  Porous  matrix  temperatures 

Figs.  3-5  compare  the  distributions  of  the  fluid-phase  tem¬ 
perature  (0f)  and  the  solid-phase  temperature  (0S)  in  the  porous 
cathode  under  various  Stanton  numbers.  The  Reynolds  number 
and  Darcy  number  are  fixed  at  Re=  12,  and  Da  =  9.8  x  10-6, 
respectively. 

A  low  fluid-phase  temperature  from  the  module  inlet  (X  =  1 .0) 
increases  as  the  flow  approaches  the  catalyst  layer  ( X = 0)  and 
then  moves  downstream.  In  contrast,  the  solid-phase  temperature 
decreases  along  X-direction  from  the  catalyst  layer  to  the  rib- 
shoulder  surfaces.  Near  the  module  inlet,  the  hot  solid  matrices 
transfer  a  part  of  heat  to  the  cold  fluids  via  convection.  Around 
the  rib  shoulder,  the  fluid-phase  temperatures  are  higher  than 
the  solid-phase  ones.  This  indicates  that  the  heat  possessed  with 
the  high-temperature  fluids  not  only  is  carried  downstream  by 
bulk  motion  but  is  also  convected  to  the  solid  matrices,  and  is 
then  conducted  out  through  the  rib- shoulder  surfaces.  When  the 
Stanton  number  increases  to  St=  14.73,  as  show  in  Fig.  4,  the 
temperature  discrepancy  between  the  solid  and  fluid  phases  is 
less  significant.  It  means  that  the  porous  medium  approaches 
the  local  thermal  equilibrium.  As  the  Stanton  number  further 
increases  to  St=  1.47  x  103,  as  shown  in  Fig.  5,  the  distribu¬ 
tions  of  the  fluid-phase  and  solid-phase  temperatures  are  almost 


Min:  1.00 


Fig.  5.  Comparison  of  fluid-phase  and  solid-phase  temperature  distributions 
inside  the  porous  electrode  for  Re  =  6,  Da-  9.8  x  10'6,  and  St=  1.47  x  103. 
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Fig.  6.  Effect  of  heat  transfer  coefficient  on  the  temperature  difference  between 
the  solid  phase  and  fluid  phase:  (a)  hy  =  5.0  x  105;  (b)  hy  =  2.0  x  106,  and  (c) 
hv  =1.0  x  109  W  m-3  K. 


the  same,  except  for  the  region  near  the  module  inlet  (Y=0.5, 
0.9  <X<  1.0).  A  large  similarity  in  the  temperature  distribution 
between  the  solid  matrices  and  the  fluids  gives  an  indication  that 
the  thermal-fluid  field  in  the  porous  electrode  has  become  local 
thermal  equilibrium. 

Fig.  6  further  shows  the  effect  of  hy  (5.0  x  105,2.0  x  106,and 
1.0  x  109  Wm~3  K-1)  on  the  temperature  difference  between 
the  solid  matrices  and  the  reactant  fluids.  Clearly,  the  solid  matri¬ 
ces  and  reactant  fluid  have  the  same  temperatures  in  the  catalyst 
layer  (— 0.25  <  A<0).  In  the  diffusion  layer  (0<X<  1.0),  how¬ 
ever,  the  temperatures  between  the  reactant  fluids  and  the  carbon 
fibers  are  different.  At  the  module  inlet,  the  fluid-phase  tem¬ 
perature  is  lower  than  that  of  the  solid-phase.  In  contrast,  the 
solid-phase  temperature  is  lower  than  the  fluid-phase  one  near 
the  rib  shoulder.  The  maximum  temperature  difference  is  about 
—01 15  foxhy  =  5.0  x  1 05  W  m  3  K  1 .  It  decreases  with  increas¬ 
ing  the  heat  transfer  coefficient. 

3.3.  Heat  fluxes 

Fig.  7  shows  the  heat  flux  distributions  along  the  X  =  0.4 
through  which  all  heat  should  be  removed  convectively  by 
fluid  (<2f)  or  conducted  out  by  solid  ((9S).  The  non-dimensional 
heat  flux  through  the  solid  and  fluid  phases  on  the  reaction 
surfaces  can  be  represented  by  Qs  =  —  (d0s/dX)w  and  Qf  = 
—  Rk(d6f/dX)w,  respectively.  It  is  seen  that  the  heat  transfer  by 
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Fig. 7.  Heat  flux  distributions  along  the  X=  1.0  for  Da  =  9.8  x  10  6,and7fe=12. 


the  solid-phase  conduction  is  higher  than  that  by  the  fluid  phase. 
This  is  because  the  significantly  higher  thermal  conductivity  of 
the  solid  phase  (Table  1)  can  provide  a  wider  thermal  pathway 
for  channelling  out  the  heat. 

3.4.  Species  fluxes 

The  oxygen  fluxes  and  the  water- vapour  fluxes  in  the  porous 
cathode  are  shown  as  vectors  in  Figs.  8  and  9,  respectively. 
The  total  species  fluxes  can  be  decomposed  into  two  parts,  i.e., 
convective  fluxes  and  diffusive  fluxes.  As  shown  in  Fig.  8(a), 
the  vectors  of  oxygen  convective  fluxes  direct  from  the  module 
entrance  toward  the  module  outlet,  which  are  similar  to  those 
of  the  velocity.  As  for  the  oxygen  diffusive  fluxes  (Fig.  8(b)), 
the  vectors  direct  from  the  module  entrance  toward  the  catalyst 
layer  and  decrease  in  magnitude  gradually.  This  is  because  the 
oxygen  is  consumed  by  the  electrochemical  reaction  in  the  cat¬ 
alyst  layer.  It  is  noted  that  the  oxygen  diffusive  fluxes  are  larger 
than  the  oxygen  convective  fluxes.  The  water- vapour  vectors 
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Fig.  8.  Comparison  of  convective  and  diffusive  species  fluxes  of  oxygen  in  the  porous  cathode  of  a  polymer  electrolyte  fuel  cell:  (a)  convective  fluxes;  (b)  diffusive 
fluxes. 
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Fig.  9.  Comparison  of  convective  and  diffusive  species  fluxes  of  water  vapour  in  the  porous  cathode  of  a  polymer  electrolyte  fuel  cell:  (a)  convective  fluxes;  (b) 
diffusive  fluxes. 


1130 


C.H.  Chao,  A.J.J.  Hwang  /  Journal  of  Power  Sources  160  (2006)  1122-1130 


of  convective  fluxes  (Fig.  9(a))  are  largely  similar  to  those  of 
the  oxygen  convective  fluxes  (Fig.  8(a)).  The  water  vapours  are 
produced  by  the  electrochemical  reaction  in  the  catalyst  layer, 
and  thus  direct  their  fluxes  diffusively  from  the  catalyst  layer 
toward  the  module  entrance.  Again,  the  water- vapour  diffusive 
fluxes  are  significantly  larger  than  the  water- vapour  convective 
fluxes. 

4.  Conclusions 

Detailed  phase-temperature  distributions  inside  a  porous 
cathode  of  a  polymer  electrolyte  fuel  cell  in  contact  with  a 
typical  interdigitated  flow  field  have  been  studied  numerically. 
The  thermal  interaction  between  the  solid  matrices  and  the 
reactant  fluid  in  the  porous  electrode  of  a  polymer  electrolyte 
fuel  cell  is  discussed  in  detail  by  examining  the  effect  of  the 
intrinsic  heat  transfer  coefficient.  A  two-equation  local  thermal 
non-equilibrium  model  is  implemented  to  investigate  the  heat 
transfer  into  the  fuel  cell  porous  cathode.  In  species  transports, 
the  Bruggemann  model  is  employed  to  describe  the  effective 
diffusivities  of  the  oxygen  and  water  vapour  in  the  porous  elec¬ 
trode. 

Results  show  that  the  wall  temperature  decreases  with 
increasing  the  intrinsic  heat  transfer  coefficient.  The  local  max¬ 
imum  wall  temperature  always  occurs  at  the  downstream  end  of 
the  module.  As  the  intrinsic  heat  transfer  coefficients  increase, 
the  porous  electrode  becomes  local  thermal  equilibrium  with  a 
strong  thermal  interaction  (heat  transfer)  between  the  solid  and 
fluid  phases.  Under  the  conditions  of  high  intrinsic  heat  transfer 
coefficients,  the  temperature  difference  between  the  solid  matri¬ 
ces  and  the  reactant  fluids  are  negligible. 
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